fid=fopen('groups.dat','r')

tline = fgets(fid);
gnum=0;

while ischar(tline)
    gnum=gnum+1;
%    area{gnum}=tline;
    areas{gnum}=tline;
    tline = fgets(fid);
    celltype{gnum}=tline(1:end-1);
    tline = fgets(fid);
    t=sscanf(tline,'%d %d %d %d');
    la(gnum)=t(1); lb(gnum)=t(2); neuron_id(gnum,1)=t(3); neuron_id(gnum,2)=t(4);
    tline=fgets(fid);
end;